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Abstract Ratios of the ip' over the J/ip production cross sections in the dilepton channel for 
C, Ti and W targets have been measured in 920 GeV proton-nucleus interactions with the 
HERA-B detector at the HERA storage ring. The ip' and J/ip states were reconstructed in 
both the and the e + e~ decay modes. The measurements covered the kinematic range 

—0.35 < xf < 0.1 with transverse momentum pr < 4.5GeV/c. The ip' to J/ip production 
ratio is almost constant in the covered xf range and shows a slow increase with pr- The 
angular dependence of the ratio has been used to measure the difference of the ip' and J/ip 
polarization. All results for the muon and electron decay channels are in good agreement: 
their ratio, averaged over all events, is R^i (n) / R^i (e) = 1.00 ± 0.08 ± 0.04. This result 
constitutes a new, direct experimental constraint on the double ratio of branching fractions, 
(B'(fi) ■ B(e)) I (B(fi) ■ B'(e)), of ip' and J/ip in the two channels. 



1 Introduction 

The hadroproduction of charmonium can be described by models [1-4] which rely on measured 
structure functions to describe partons in the initial state, perturbative QCD to calculate the 
production of intermediate cc states, and a subsequent non-perturbative hadronization step. 
Depending on the model assumptions and scope, this latter step must account for several effects, 
including color neutralization by soft gluon emission, bound state formation and the influence of 
nuclear matter ( "nuclear suppression" ) . This in turn leads to the introduction of free parameters 
which are tuned to describe existing data. The tuned model can then be used to predict the 
production of charmonium states not considered in the tuning step or in other kinematic regimes. 

Production of the ip' in proton induced collisions has been measured over a wide range of 
energies in several experiments at Fermilab and CERN. But the experimental situation is still not 
satisfactory since most measurements suffer either from limited sample size or from additional 
large uncertainties resulting from signal extraction in the presence of large backgrounds. In many 
of the fixed-target experiments, absorber material was placed behind the target to reduce the 
muon background produced by pion and kaon decays in flight. For such experiments the e + e~ 
mode was not observable and also the mass resolution suffered from the scattering of the muons 
in the absorber. In some cases, the mass resolution was insufficient for a separation of the J/ip 
and ip' signals. 

In contrast, using the HERA-B detector's advanced particle identification and tracking sys- 
tems, it was possible to reduce the background from ir and K decays in flight by using track 
quality criteria without inserting absorber material, and thereby to have good mass resolution 
and also afford the possibility of triggering on and reconstructing electrons. The pi + ji~ and e + e~ 
decay channels were recorded simultaneously. The rather large HERA-B samples of J/ip and ip' 
(300,000 and 5000 events , respectively) and well separated J/ip and ip' signals will contribute 
to the clarification of the experimental situation. 

In order to minimize possible systematic biases introduced by acceptance and efficiency cor- 
rections, we report the ratio of the ip' and J/ip production cross section times the dilepton 
branching ratio. The sample size allows measurements of this ratio in several bins of the Feyn- 
man variable, xf, and transverse momentum, px- The decay angular distribution of the ratio, 
which is related to the difference of the ip' and J/ip polarization was also studied. 

The branching ratio of ip' in the dimuon decay channel is not well known. The PDG value 
of the branching ratio B'(/j,) = (0.73 ± 0.08)% [5] for ip' has a relative error of 11% and is 
a factor of 3 less accurate than the corresponding branching ratio for the e + e~ decay, B'(e). 
Furthermore almost all measurements of ip' to J/ip production ratio have been performed in the 
pt + pt~ channel only. The production ratio measured in pp collisions by detecting e + e" pairs [6] is 
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rather imprecise (0.019 ± 0.007). Our relatively large samples in both channels therefore permit 
a sensitive test of e — /i, universality in ip' dilepton decays. 

The paper is structured as follows. After a short description of the HERA-B detector, the 
trigger and the data sample in Sect. 2, the method used for the measurement of the ip' and 
J/ip production ratio is introduced in Sect. 3. The Monte Carlo simulation and the event re- 
construction and selection are presented in Sects. 4 and 5. In Sect. 6, the difference of ip' and 
J ftp polarization is studied, followed by a discussion of sources of systematic uncertainties in 
Sect. 7. Finally, the experimental results are presented, the differential ijj' and J/ip production 
ratio in Sect. 8, the inclusive ip' and J/ip production ratios in Sect. 9 and the nuclear dependence 
in Sect. 10. Sect. 11 confronts the HERA-B results with previous measurements and theoretical 
model predictions. The paper ends with conclusions given in Sect. 12. 

2 Detector, trigger and data sample 

The fixed target experiment HERA-B [7,8] was operated at the HERA storage ring at DESY 
until 2003. The forward spectrometer of HERA-B covered angles from 15 to 220 mrad in the 
bending plane and from 15 to 160 mrad vertically. The target [9] consisted of two independent 
stations, separated longitudinally by 4 cm along the beam. Each station contained four wires 
which could be independently positioned in the halo of the HERA proton beam (inner, outer, 
below and above with respect to the beam). Various materials (C, Ti, W) were used for the 
wires. 

The tracking system consisted of the Vertex Detector System (VDS) followed by a magnetic 
dipole field of 2.13 Tm and a main tracker. The VDS [10] was a silicon strip detector with a 
pitch of 50 /im. The Inner Tracker (ITR) [11] of the main tracking system used micro-strip gas 
chambers with 300 /im pitch in the inner part of the main tracker up to distances of 25 cm around 
the proton beam. The Outer Tracker (OTR) [12,13] used honeycomb drift chambers with 5 and 
10 mm diameter cells covering the remaining acceptance. A spatial resolution of about 150 /im 
and 360 /im was achieved for the ITR and OTR, respectively. 

The main tracker system consisted of four stations immediately behind the magnet: the 
Pattern Chambers (PC), and two stations positioned further downstream: the Trigger Chambers 
(TC). A Ring Imaging Cherenkov detector (RICH) [14] was situated between the PC and TC 
stations. The TC chambers were mainly used by the trigger. The TC also served together with 
the PC to provide additional measurements of the track position and direction. The number of 
layers was 30 and 12 in the PC and TC chambers, respectively. 

Particle identification was provided by the RICH, a MUON system and an 
Electromagnetic Calorimeter (ECAL). The ECAL [15] followed by the MUON system were 
located at the end of the spectrometer. The ECAL was optimized for good electron/gamma 
energy resolution and electron-hadron discrimination. The ECAL had variable granularity: the 
cell sizes were 2.2 cm, 5.5 cm and 11cm in the inner, middle and outer part, respectively. The 
ECAL was instrumented with a fast digital read-out and a pretrigger system. The MUON 
system [16] consisted of four tracking stations at various depths in an iron or concrete absorber. 
Only muons above 5 GeV/c have a significant probability of penetrating through the absorbers. 

The trigger chain included the ECAL and MUON pretrigger systems which provided lepton 
candidate seeds for the First Level Trigger (FLT) which in turn required at least two pretrigger 
seeds. Starting from the seeds, the FLT attempted to find tracks in a subset of the OTR tracking 
layers, requiring that at least one of the seeds results in a track. Starting once again from the 
pretrigger seeds, the Second Level Trigger (SLT) searched for tracks using all OTR layers and 
continued the tracking through the VDS, finally requiring at least two fully reconstructed tracks 
which were consistent with a common vertex hypothesis. The pretriggers and FLT are hardware 
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triggers which reduce with an output rate of typically 25 kHz while the SLT is a software trigger 
running on a farm of 240 Linux PCs. 

This analysis is based on 164 million dilepton triggered events collected in the 2002-2003 
physics runs. About 300,000 J/ip lepton decays contained in the data sample are almost equally 
divided between the dimuon and dielectron channels. The wire materials used are carbon (A = 
12, about 65% of the J/ip data), tungsten (A = 184, about 31%) and titanium (A = 48, about 
4%). Only runs in which detector components performed well and trigger conditions were stable 
were used. Six periods with constant experimental conditions are defined and the data are 
grouped accordingly. 



3 Measurement method 

The analysis is based on the selection of dilepton events and fitting of the dilepton invariant 
mass spectra in the area around the J/ip and ip' signals. The number of events in the J/ip peak, 
N^p, is given by: 

Nj, = a(J/ip) ■ B(J/ip -» l + l~) ■ C ■ e , 

i.e. the product of the cross section (a(J/ip)), the branching ratio into dilepton pairs (B), 
the integrated luminosity (£) and the total reconstruction efficiency (e). The efficiency takes the 
following effects into account: the probability for leptons to be within the detector acceptance and 
to be properly registered in the detectors, the probability for the dilepton pair to be triggered, the 
probabilities to reconstruct the tracks and the dilepton vertex, and the probability to select the 
dileptons for the final analysis. The luminosity is identical for all charmonium states produced 
on the same targets, and therefore cancels in the ratios. The ratio of the ip' and J/ip cross 
sections in the l + l~ channel, R^'(l), is equal to: 

where I denotes the leptonic decay channel (e or /i), a(a') is the J/ip (ip 1 ) production cross 
section and B (B 1 ) is the branching ratio for the l + l~ decay of the J/ip (ip') meson. 

The ratio N^i / or "raw ip' /ip ratio", is defined from the fit of the J/ip and ip' signals and 
must be corrected by the efficiency ratio, e/e', where e is defined above and e' is the corresponding 
efficiency for detection of ip' mesons. The efficiencies are evaluated by a Monte Carlo simulation. 

For the analysis we select the kinematic domain 

-0.35 < x F < 0.10 and < p T < 4.5 GeV/c (2) 

for the Feynman variable xf and the transverse momentum pt, respectively. This domain cor- 
responds to our acceptance range in xf and at high px is limited by lack of ip' events, i.e. is 
determined by the sample size. 



4 Monte Carlo simulation 

Results of Monte Carlo simulations (MC) were used to determine the efficiencies of J/ip and 
ip' triggering, reconstruction and selection. J/ip and ip' production was simulated by generating 
the basic process, pN — > ccX with PYTHIA 5.7 [17] and hadronizing the c and c quarks with 
JETSET 7.4 [17]. The momentum and energy of the remaining hadronic system, X, is given 
as input to FRITIOF 7.02 [18] which generates particles in the underlying event taking into 
account their interactions inside the nucleus. 
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The generated ip' events were assigned weights such that the resulting weighted xf distribu- 
tion of the if)' matched that of the J/ip, since both measurements and theoretical models show 
little or no difference of these distributions in the kinematic domain Eq. (2) as described in 
Sect. 8. Additional weights were applied to both J/ip and ip' events such that the resulting pr 
distribution matches 

-6 

, (p T ) = 1.29GeV/c. (3) 

The average transverse momentum, (pt), m (3) was taken from a preliminary HERA-B analysis 
of the J/ip sample [19]. The effects of possible differences in pt and xp distributions for the J/ip 
and ip' are described in Sect. 7. 

A GEANT 3.21 [20] based package [21] performed tracking of the particles through the 
HERA-B detector and the simulation of the detector response. The status (dead/alive), efficiency 
and noise level of each detector channel were inferred from the data separately for each of six 
running periods (see Sect. 2) and used in modeling the detector response. The MC events were 
reconstructed with the same code used for the data. 

5 Event reconstruction and selection 

The track reconstruction included the following steps: finding straight track segments in the 
VDS and PC area, propagation of the PC segments through the TC area, matching of VDS 
and PC segments and a full iterative fit of the tracks. The track segments were propagated into 
the TC area to provide a better measurement of the tracks in the RICH, ECAL and MUON 
detectors. Particle identification estimators were evaluated for each track using information from 
these three detectors. 

Pairs of like sign and unlike sign candidates of muons and electrons were selected. A vertex 
fit was performed for each candidate with a weak cut on the x 2 probability (> 10 -5 ) to ensure 
that the tracks originated from the same vertex. At this stage, 99% of the recorded J/ip and if)' 
candidates with reconstructed lepton tracks were selected. Further cuts and fits of the dilepton 
invariant mass spectra which were specific to the muon and electron modes will be discussed in 
the following two subsections. 

5.1 Selection of J/ip, if)' — > 

For the selection of fi + /J,~ pairs, soft cuts were applied to the momenta (400GeV/c > p > 
6GeV/c) and transverse momenta (px > 0.7GeV/c) of the muons. The lower cut values were 
close to the trigger requirements and the upper momentum cut rejected tracks outside the 
expected range for J/ip(ip') decay products. 

The main background in the muon channel was from muons produced in K, tt — > [iis decays 
in flight. If a decay occurred inside the tracking system, the "broken trajectory" could have led 
to a large x 2 i n t ne track fit. A decay between the VDS and the PC could additionally have led 
to a poor match between VDS and PC track segments. The tracks from the main tracker were 
projected into the MUON system, where hits were associated with the track and used to obtain 
a x 2_ value with respect to the extrapolated track, taking into account multiple scattering. This 
muon-x 2 was used to obtain a quality estimator varying from to 1 and called the "muon 
likelihood". Kaon and pion decays between the PC and the MUON system could have led to 
a large muon~x 2 and therefore a low muon likelihood, as could wrongly assigning MUON hits 
produced by another particle. 



da 
dp\ 
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Figure 1. The invariant mass distribution of n + fi~ in the region of the J/ip and ip' peaks for selected dimuons 
from all three targets and fitted with the sum of functions describing the signals and the background (dashed 
curve) . 

Thus, muon candidate tracks were selected according to the x 2 °f their track fits, qualities 
of their VDS-PC matches and their muon likelihoods in such a way that only about 5% of 
real muons were discarded. Tracks were also required to have a low probability of being a kaon 
according to the RICH analysis. 

With these selections, the background under the J/ip and ip' peaks was reduced by a factor 
of 2.5. The signal loss due to these selections was distributed rather evenly over the kinematic 
range defined by Eq. (2) and was reproduced by the MC. After the selection, (89.4 ± 0.2)% of 
the events in the J/ip peak survived in the data and (89.1 ± 0.1)% in the MC. 

The dimuon mass spectra for both J/ip and ip' were described as a sum of two functions [22]: 
1) a symmetric function, being a superposition of three Gaussians, which takes into account 
track resolution and effects of Moliere scattering, 2) a function representing a radiative tail due 
to the emission of a photon in the final state of the dimuon decay J/ip(ip') — > /j, + fi~'y. 

The presence of the radiative tails and tails from Moliere scattering in the invariant mass 
distribution of the J/ip and ip' signals demanded that both peaks be fit simultaneously. The total 
number of events in a peak was obtained by an integration of the fitted function over the full 
range of possible dimuon masses. The mass of the ip 1 was obtained by scaling the fitted J/ip mass 
by the factor Mo(ip') /Mo (J ftp) , where the Mq were the nominal values of the respective masses. 
The fitted J ftp mass resolution was rescaled for the tp' using the dependence of the momentum 
resolution for the HERA-B spectrometer. The background was described as an exponent of a 
quadratic polynomial of the dimuon mass. 

The functions describe the dimuon mass spectra well as shown in Fig.l, where x 2 = 154 
for 153 degrees of freedom. The fitted mass of J/ip is 3093 MeV/c 2 and the mass resolution 
(FWHM / 2.35) is 38 MeV/c 2 . 

5.2 Selection of J/ip, ip' — ► e+e - 

In the dielectron channel, the background was larger than for muons mainly due to misidenti- 
fication of charged pions interacting in the ECAL and overlaps between energy deposits from 
photons and charged hadrons. Despite the higher background a careful study of the electron 
identification cuts resulted in J/ip and ip' signals of comparable purity and significance. 
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The first step of the event selection consisted of the requirement that track momenta lie 
between 4 and 400 GeV/c, corresponding to the expected range for J/ijj and ip' decay products. 
Finally, an additional cut on the transverse energy of the ECAL cluster {Et > 1.15 GeV) 
was applied in order to equalize different cut thresholds used at the pretrigger level during 
the acquisition periods. The signal selection was improved by rejecting electron/positron pairs 
having a distance of closest approach between the two tracks greater than 320 //m. 

The reconstructed momentum vectors of electrons and positrons were corrected for energy lost 
due to bremsstrahlung (BR) emission in the material before the magnet (evaluated as about 
18% on average). BR was identified and determined for each track by looking for an energy 
deposition in coincidence with the extrapolation of the track's VDS segment to the ECAL. The 
recovered energy was then added to the momentum measured by the tracking system. Being a 
clear signature of electrons, BR emission was also exploited to obtain a substantial background 
reduction, essential for an accurate counting of the small ip' signal. It was therefore required 
that at least one lepton of the decaying J/tp or ip' had an associated BR cluster (defined as an 
ECAL cluster close enough to the VDS segment extrapolation). 




Figure 2. The invariant mass distribution of e + e~ in the region of the J/tp and tp' peaks for selected dielectrons 
from all three materials, fitted with the sum of the functions describing the signals and the background (dashed 
curve) . 



Two additional types of selection cuts were applied to those electrons and positrons which 
had no associated BR cluster. The first of these cuts is on the E/p ratio, where E is the energy of 
the electron/positron cluster and p is the associated track momentum. The E/p distribution for 
electrons has a Gaussian shape with mean value close to 1 and width, / p , varying between 6.4% 
and 7.4% depending on calorimeter section. Values of E/p far from 1 correspond to particles, 
mostly hadrons, which release only part of their energy in the calorimeter. The second cut was 
applied to the distance between the reconstructed cluster and the track position extrapolated to 
the ECAL {A x and A y for the x and y direction). The A x and A y distributions for electrons are, 
apart from a small tail, well described by Gaussians centered at zero, with widths between 0.2 cm 
and 1.0 cm, depending on calorimeter section. A cut on these quantities leads to a significant 
reduction of the contamination from hadrons and random cluster-track matches, for which the 
^x(y) distributions are expected to be significantly wider compared to those of the electrons. 
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All the above requirements were simultaneously optimized by maximizing the significance of 
the ip' MC signal. The accepted ranges for E/p, A x and A y were respectively determined as 
— 1.0 ■ <Je/p < E/p — 1 < 3.2 • (Je/p, \Ax\ < 1.7 a x and \ Ay\ < 1.8 a y . As a result, the S/B ratio, 
evaluated for the J/ ip (the only one of the two charmonium states visible also before the selection 
cuts), increased by about a factor of 10 with respect to the triggered events, with an efficiency 
for signal selection of (34 ± 3)% in the real data and (35.1 ± 0.4)% in the MC simulation. 

To count the number of J /ip and ip', a Gaussian for the right part of the peak and a Breit- 
Wigner for the left part were used, allowing for a sizeable asymmetry of the signal due to the 
not fully reconstructed BR emission and a contribution from J/ip — > e + e~7 decay. Both the 
width and the asymmetry of the ip' signal were kept at a fixed ratio (determined from the MC) 
with respect to those of the J/ip. The background was parametrized as a Gaussian for the low 
invariant mass values and an exponential at higher mass, requiring continuity of the resulting 
function and its first derivative. The invariant mass distribution plotted in Fig. 2 shows our final 
selection of 68500±470 J/ip and 1680±130 ip': the resulting mass and width (of the Gaussian) of 
the J/ip are, respectively, 3110 and 72 MeV/c 2 . The \ 2 of the fit is 89 for 69 degrees of freedom. 



6 Angular dependence of R^> 

Measurements of polarization (the commonly used term to denote spin alignment) provide one 
of the most significant tests of models of charmonia production [23]. Understanding polarization 
at moderate and low pt is also important since most of the published ^-integrated data come 
from the lower p? region. The polarization of charmonium is measured by observation of the 
angular distribution in its decay into 

^ oc (1 + (4) 

where 9 is the polar angle of the / + with respect to the z axis of a polarization frame defined in 
the rest system of the charmonium state. In general, the parameter A depends on the definition 
of the polarization frame [23]. The commonly used polarization frames are specified by the choice 
of the z axis. The three-momentum of the beam particle, Pb, and that of the target, pt, in the 
rest system of the charmonium can be used to define the z axis. 

In the Gottfried-Jackson frame, the z axis is parallel to pb, in the s-channel helicity (or recoil) 
frame, the z axis is defined as the direction of the charmonium three-momentum in the hadronic 
center-of-mass frame, i.e. it is parallel to —(pb + Pt), and in the Collins-Soper frame, the z axis 
bisects the angle between pb and —pt- 

For J/ip production, several measurements of polarization have been published. The parame- 
ter A measured in pBe collisions [24] in the Gottfried-Jackson frame is A = 0.01 ±0.12 ±0.09 and 
A = -0.11±0.12±0.09 at 530 GeV and 800 GeV, respectively. A high statistics measurement [25, 
26] in 800 GeV pCu interactions (Collins-Soper frame) resulted in A = 0.069 ± 0.004 ± 0.08. The 
quoted errors are statistical and systematic. A measurement [27] for ip' was performed for ir~W 
interactions at 253 GeV where A' = 0.02 ± 0.14 was obtained in the Gottfried-Jackson frame. 
All results are consistent with a small difference in the ip' and J/ip polarization. The collider 
result [28] (s-channel helicity frame) A' = —0.08 ± 0.63 ± 0.02 for a mean p? equal to 6.2 GeV/c 
is also consistent with that, but statistically less significant and, because of the relatively high 
Pt, not directly comparable with the fixed-target results. 

In the kinematic domain (2) of HERA-B, the ip' signal is clearly visible in the limited interval 
| cos 9\ < 0.6(0.8) for electrons (muons) in the Gottfried-Jackson frame. In the s-channel helicity 
frame the acceptance is more uniform as a function of cos 6 and R^i is measurable over the full 



10 



0.03- 



■o 
ta 

"e 

T3 

CD 0.02 



0.01 




cose 



Figure 3. i?^/ as a function of cosB in the s-channel helicity frame. The dashed curve is the fitted function (5) 
for A = and the fitted parameter A' = 0.23 ± 0.17. 



Table 1. The ratio of tp' and J/tp cross sections in the muon, R^i(fi), electron, R^>(e), channels and for both 
channels combined, R^./ , obtained for different cos 8 intervals using the data from all targets. The value 9 is the 
polar angle of the positive lepton in the s-channel helicity frame. Additional systematic uncertainties of 5% for 
the electron channel and 3% for both the muon channel and for the combined result (see Table 2) are not included 
here. 



cos 8 


R^' {n) 




R-ip' 


-1.0- 


--0.8 


0.0148 ± 0.0024 


0.0144 ± 0.0054 


0.0147 ±0.0022 


-0.8- 


--0.6 


0.0211 ± 0.0022 


0.0169 ± 0.0046 


0.0203 ± 0.0020 


-0.6- 


--0.4 


0.0203 ± 0.0019 


0.0107 ± 0.0034 


0.0180 ±0.0017 


-0.4- 


--0.2 


0.0154 ± 0.0016 


0.0185 ± 0.0033 


0.0160 ± 0.0014 


-0.2- 


- 0.0 


0.0152 ±0.0016 


0.0130 ± 0.0028 


0.0147 ±0.0014 


0.0- 


- 0.2 


0.0161 ±0.0016 


0.0140 ± 0.0028 


0.0156 ±0.0014 


0.2- 


- 0.4 


0.0152 ±0.0017 


0.0159 ±0.0030 


0.0154 ±0.0015 


0.4- 


- 0.6 


0.0187 ±0.0019 


0.0167 ±0.0035 


0.0182 ±0.0017 


0.6- 


- 0.8 


0.0157 ±0.0021 


0.0154 ± 0.0043 


0.0156 ±0.0019 


0.8- 


- 1.0 


0.0176 ±0.0026 


0.0269 ± 0.0060 


0.0191 ± 0.0024 



interval of cos 9. We therefore prefer to present our results in the s-channel helicity frame. The 
values of R^i for ten cos# bins are shown in Table 1. 
The dependence of R^i on cos 6 can be described as 

B 1 da' /d cos 6 1 + A'cos 2 # 

B da Id cos 9 * 1 + Acos 2 ^' ^ ' 

where A' and A are related to the polarization of the ip' and the J/i/j, respectively. The fit of 
R^> (see Fig. 3) with this function, where A' is a free parameter and A fixed, depends on the A 
value, but the difference A' — A is nearly constant. We obtain in the s-channel helicity frame 

AX = A'- A = 0.23 ±0.17, (6) 

for the fit with A = 0. To check the stability of the result, we varied A in the interval (—0.2, 
0.2), which is broader than the limits given by the published measurements of J /ip polarization, 
and obtained a variation of AX in the range ±0.04, which is much smaller than the statistical 
uncertainty. 
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7 Systematic errors of R^> 

The correction of the raw ip' /ip ratio in (1) using the efficiency ratio, e/e', calculated by MC was 
the subject of a detailed study. The efficiency ratio, e/e', was mainly determined by geometric 
factors since effects such as detector and trigger efficiencies nearly cancel out in the ratio. The 
geometrical factors can be calculated reliably with MC and were found to be quite stable. This 
is confirmed by the data since the raw ip' /ip ratio is also quite stable for different periods of data 
taking and target configurations. 

Nevertheless, a small systematic bias can appear due to uncertainties in the simulation of 
kinematical distributions of charmonia production as well as of the detector and trigger. Possible 
biases were evaluated by varying parameters of the MC simulation. 

The MC simulation assumed unpolarized production for J/ip and ip'. However, due to the 
limited acceptance, the efficiency ratio e/e' in Eq. (1) could depend on a difference between the 
simulated decay angular distributions for the ip' and J/ip. The ratio e/e' depends slightly, almost 
linearly on the difference A' — A and consequently 

4^ = 0.2 • AX. (7) 

To obtain Eq. 7, we evaluated the efficiency ratio e/e' by simulating ip' and J/ip production for 
different AX values in the Gottfried-Jackson frame. In this frame, the acceptance is less uniform 
as a function of cos 6, and the bias of R^i due to the possible difference of angular distributions 
for the ip' and J /if; is dominated by the size of AX. Theoretical arguments suggest that the 
difference in the ip' and J/ip polarization is small. The color evaporation model (CEM) predicts 
that A and A' are equal [29]. In the Nonrelativistic QCD (NRQCD) approach, the polarization of 
the J/ip is larger than that of the ip' due to feeddown from the % c , but A is unlikely to exceed A' 
by more than 0.2 [23]. For the HERA-B experiment, the values for A calculated in the NRQCD 
approach [29] range from to 0.1. We estimated the contribution to the relative systematic error 
due to polarization a sys {R^i ) /R^i = 0.02, corresponding to AX = 0.1 in (7). 

Theoretical predictions for J/ip and ip' polarization in fixed-target experiments are largely 
independent of the target and beam types and vary slowly with beam energy [29]. Therefore 
direct comparison of results of different experiments is meaningful. Combining measurements 
[24,25,30] for A, [27] for A' and the HERA-B result (6), we obtained AX = 0.09 ± 0.11. The 
latter estimate disfavors also a big difference between A and A'. 



Table 2. Estimated relative systematic errors of the R^i in the n + fi , R^r(n), and e + e , R^,i(e), channels from 
various sources, as listed. The total uncertainty was obtained by quadratic summation. 



Source of uncertainty cr sys (J?^')/-^V'' (m) o"s ys (-R^') / 'R^i(e) 



Polarization 


0.02 


0.02 


Pt distribution 


0.0006 


0.0006 


xf distribution 


0.004 


0.004 


Trigger simulation 


0.006 


0.02 


Counting method 


0.014 


0.035 


Total 


0.03 


0.05 



To estimate the importance of the MC model used for the production kinematics, the param- 
eters of the generated xp and pr distributions were changed within variations allowed by the 
measured differences between ip' and J ftp kinematics [19,31,32]. The corresponding variations of 
R^i were found to be rather small as listed in Table 2. Contributions to the systematic error due 
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to differences of kinematical distributions and polarization of tp' and J/ ip were highly correlated 
for muons and electrons. 

Uncertainties connected to the trigger were evaluated by removing the trigger simulation for 
both muons and electrons and, as a second test, increasing the value of the Et threshold with 
respect to the one used by the pretrigger for electrons. These are quite extreme tests which surely 
overestimate possible trigger biases, nonetheless, as seen in Table 2, the changes in the R^i do 
not exceed 1% and 2%, for muons and electrons, respectively. Moreover, it has been verified for 
electrons that any systematic bias due to the BR request on the value of R^i is negligible. 

The last item listed in Table 2 is related to different methods of counting the numbers of J ftp 
and tp' in (1). We obtained these values either by taking all events in the kinematic domain (2) 
or by summing the events in the bins of the xp, px or cos 8 distributions. In doing so, the 
efficiencies were evaluated either for the whole kinematic domain or for the different bins. 

The total relative systematic error on R^i was estimated to be 3% and 5% for the muon and 
electron modes, respectively. The larger error for the electron mode can be attributed to the 
higher contribution of background. 



8 x F and p T dependence of R^> 

To obtain the ratio of the tp' to J/ip production cross sections times the dilepton branching ratio, 
FL)i , as a function of xf we selected dilepton pairs in eight x f intervals and performed the full 
analysis for each interval independently. The results are listed in Table 3, separately for 
e + e~ channels as well as for the combined result. Fig. 4 displays the HERA-B results together 




-0.2 



0.2 



0.4 



0.6 



Figure 4. Measurements of R^i as a function of xf- Results of this work, combined for fi and e decay modes 
and for all targets, are shown together with previous experiments which measured charmonia production in pN 
(NA50 [33], E771 [31], E789 [34]) and 7r~N interactions [27]. Only statistical errors are shown. Two calculations, 
one using the NRQCD approach [36] and the other using the CEM model [35], are indicated by the dashed 
and dotted curves, respectively. Both calculations were made for pp interactions and corrected for the nuclear 
suppression in the HERA-B target as described in the text. 



with earlier measurements. Neither the HERA-B nor the NA50 data [33] exhibit significant xp 
dependencies. R^i measurements of E771 [31] in interactions extend up to Xi? = 0.2, but 
suffer from large statistical errors. A statistically more accurate measurement, but performed 
in 7r~N interactions [27], also indicates only a slow variation of in the range of low and 
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Table 3. The ratio of ip 1 and J/tp cross sections in the muon, electron, _R^/(e), channels and for both 

channels combined, R^i, obtained for intervals in if using the data from all targets. Additional systematic 
uncertainties of 5% for the electron channel and 3% for both the muon channel and for the combined result (see 
Table 2) are not included here. 



X F 



R^> (e) 



E 



if}' 



-0.35 
-0.25 
-0.20 



-0.15- -0.10 



-0.10 
-0.05 
0.00 
0.05 



-0.25 
-0.20 
-0.15 



-0.05 
0.00 
0.05 
0.10 



0.0309 ± 0.0076 
0.0172 ± 0.0038 
0.0156 ± 0.0022 
0.0165 ± 0.0014 
0.0162 ± 0.0011 
0.0160 ±0.0011 
0.0162 ± 0.0013 
0.0095 ± 0.0025 



0.0265 ± 
0.0244 ± 
0.0202 ± 
0.0175 ± 
0.0179 ±0 
0.0152 ±0 
0.0147 ±0 
0.0165 ±0 



0088 
0061 
0040 
0030 
0025 
0022 
0025 
0037 



0.0290 ± 
0.0192 ±0 
0.0167 ±0 
0.0167 ±0 
0.0165 ±0 
0.0158 ±0 
0.0159 ±0 
0.0117 ±0 



0058 
0032 
0019 
0013 
0010 
0010 
0012 
0021 



moderate xf values. We may expect a similarity of the xf behavior of the R^> in and 7r~N 
interactions, because of the dominance of gluon fusion in the charmonium production process. 

In the CEM [1] the differential and integrated charmonium production rates for different 
charmonium states should be proportional to each other and independent of projectile, target 
and energy. The results of the CEM for pp interactions [35] show a flat x f behavior and agree with 
our results (Fig. 4). The NRQCD calculation [36] of tp' and J /ip production in pp interactions 
at 920 GeV shows a slow variation of R^i versus x p and is consistent with our measurements as 
well. The results of the NRQCD and CEM models were both corrected for nuclear suppression 
as will be discussed in Sect. 11. 

The results for R^i for eight bins in transverse momentum are given in Table 4. R^> shows 
a tendency to increase as a function of px as seen in Fig. 5. The fit of with the ratio of 
functions of the form given in Eq. (3), where the average transverse momentum, (pt)', for ijj' is a 
free parameter and (pt) for J/tp is fixed to the value (3), gives (pt)' — (pt) = (0.08±0.03) GeV /c. 
These results are consistent with those of E771 [31] and E789 [34], although the large statistical 
errors of these two previous measurements do not allow for a decisive test. 
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Figure 5. Measurements of R^i as a function of the transverse momentum. Combined results for e and fi modes 
obtained for all three targets are shown together with previous results from experiments E771 [31] and E789 [34]. 
Only statistical errors are shown. A fit of the HERA-B results with the ratio of functions (3) as described in the 
text, is shown by the dashed line. 
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Table 4. The ratio of ip' and J/ip cross sections in the muon, _R^,/(/i), electron, R^,/(e), channels and for both 
channels combined, R^i , obtained for intervals in pr using the data from all targets. Additional systematic 
uncertainties of 5% for the electron channel and 3% for both the muon channel and for the combined result (see 
Table 2) are not included here. 



p T [GeV/c] R^'{n) R^'(e) R^,, 



0.0^ 


0.5 


0.0143 ± 0.0015 


0.0129 ± 0.0031 


0.0140 ±0.0014 


0.5^ 


0.9 


0.0156 ± 0.0012 


0.0145 ± 0.0022 


0.0153 ±0.0011 


0.9^ 


1.3 


0.0167 ±0.0012 


0.0159 ± 0.0025 


0.0166 ±0.0011 


1.3-^ 


1.7 


0.0171 ± 0.0014 


0.0201 ± 0.0027 


0.0177 ±0.0012 


1.7-=- 


2.1 


0.0156 ± 0.0019 


0.0134 ± 0.0034 


0.0151 ± 0.0017 


2.1 + 


2.5 


0.0192 ± 0.0027 


0.0145 ± 0.0040 


0.0177 ± 0.0022 


2.5-=- 


2.9 


0.0248 ± 0.0041 


0.0150 ± 0.0055 


0.0213 ± 0.0033 


2.9^ 


4.5 


0.0277 ± 0.0043 


0.0198 ± 0.0061 


0.0251 ± 0.0035 



9 Combined results for R^> and for the double ratio of leptonic branching fractions 

An independent analysis was performed for each target sample (C, Ti, W) and separately for 
both dimuon and dielectron trigger modes. For each sample the dilepton invariant mass spectra 
were fitted to obtain the raw ip'/ij) ratio which was then corrected according to (1) by the 
efficiency ratio e/e'. 

The values of R^> in both dilepton channels obtained for all three targets are shown in 
Table 5. For each target the values of R^i for the muon and electron mode were combined and 
the ratio R^i (/j,) / R^r (e) was evaluated. 



Table 5. The ratio of ip' and J/ip production cross sections per nucleus in the muon, 7?^/(/i), electron, R^,(e), 
channels and for both channels combined, R^i, for each target. The ratio R^i (fx) / R^i (e) is also presented. The 
quoted errors indicate the statistical and systematic errors, respectively. 



Target 




R^'{e) 


Rip' 


R^t^/R^ie) 


C 
Ti 
W 


0.0166 ± 0.0007 ±0.0004 
0.0198 ± 0.0029 ± 0.0005 
0.0158 ±0.0011 ±0.0004 


0.0154 ±0.0012 ± 0.0007 
0.0200 ± 0.0045 ± 0.0009 
0.0182 ±0.0025 ±0.0012 


0.0163 ± 0.0006 ± 0.0005 
0.0199 ± 0.0024 ± 0.0006 
0.0162 ± 0.0010 ± 0.0005 


1.08 ±0.10 ±0.04 
0.99 ±0.27 ±0.04 
0.87 ±0.13 ±0.03 



We obtained the ratio R^* (pi) / R^' (e) , combined for all targets, as a weighted average of the 
three corresponding values from Table 5: 

iV (n) /Rip, (e) = 1.00 ± 0.08 ± 0.04. (8) 

The quoted errors are statistical and systematic. The systematic error was estimated by com- 
paring the results of different counting methods and by taking into account the uncertainty of 
trigger simulation (mainly due to the electron mode) as discussed in Sect. 7. Contributions to 
the systematic error due to differences of kinematical distributions and polarization of ip' and 
J /ip almost cancel out in (8). 

The ratio (8) also represents a direct measurement for the double ratio of the branching 
fractions of the ip' and J/i/j dilepton decays: 

= B\y) ■ B(e) 

iV(e) B<jn)-B'(ey [) 

The PDG values of branching ratios are [5]: B(ji) = (5.88 ± 0.10)%, B(e) = (5.93 ± 0.10)% for 
J/V dilepton decays and B'(ji) = (0.73 ± 0.08)%, B'(e) = (0.741 ± 0.028)% for V' decays. The 
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Figure 6. The error ellipse for B'(e) versus B\n), evaluated by using the covariance matrix from [5] (dashed 
curve) and after (solid curve) applying the HERA-B constraint given by (8). 



branching ratios are not independent since they are obtained by a global fit performed by the 
PDG [5], which results in B' (fi) / B 1 (e) = 0.99 ± 0.13. The latter value and the J/ip branching 
ratios correspond to the following determination of the double ratio: 

B'(n) ■ B(e) , , 

whM = L00±0 ' 13 ' (10) 

which is compatible with (8), but slightly less precise. 

Our result (8) for R^> (fi)/R^/ (e) would give a significant contribution as an additional con- 
straint of the global fit of ip' branching ratios as shown in Fig. 6. The error ellipse for B'(e) 
versus B'((i), evaluated by using the covariance matrix from [5], is shrunk along the B'(fi) axis 
(partially at the expense of a tilt, i.e. a correlation) after applying the additional HERA-B 
constraint by (8). 



10 The x F dependence of nuclear suppression 

Usually the nuclear dependence of cross sections is parameterized as a power law with exponent 



o 



oa oc A a , (11) 



where A is the atomic mass number. The nuclear dependence of R^/ is then also a power law in 
A: 

oc A a ' ~ a (12) 

Boa 

with the parameters a' and a for the ip 1 and J/tp, respectively. Recent high statistics measure- 
ments of the nuclear dependence of J/ip and tp' production for proton- nucleus collisions were 
made at Fermilab E866 [37] and at CERN NA50 [33]. Both results are presented in terms of 

A a = a' — a 
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Figure 7. The difference of a for ip' and J/ip production as a function of if determined from C and W data 
(circles), compared to results of E866 [37] (solid triangles) and NA50 [33] (triangles). A fit of the E866 results by 
a constant (13) is shown by the dashed line. The results of CEM for color singlet nuclear absorption [38] at 920 
GeV are displayed by the solid line. The NRQCD results [38] at 450 GeV and 920 GeV are shown by the dotted 
and dashed-dotted curves, respectively. 



for various xp slices in Fig 7. The two results are in reasonable agreement. The HERA-B 
results for A a, derived from C and W data and also shown in Fig 7, are consistent with both 
measurements, but less accurate. The measurements shown in Fig. 7 indicate little or no xp 
dependence of A a. 

The calculations of nuclear suppression by the CEM for color singlet absorption and NRQCD 
models [38] also indicate slow x p and collision energy dependences of A a and that a is larger 
than a', however the predicted difference between a' and a is less than the measured difference. 
The CEM, with the additional assumptions in [39], predicts that A a should be zero over the 
xf range of this measurement. Because of the small variation of A a in our kinematic range (2), 
we will use the experimental value averaged over xp (see Sect. 11) for the analysis of our results 
and for comparison with calculations performed for pp interactions. 

The E866 results, in the xp range shown in Fig. 7, were fitted by a constant value: 

Z\a(E866) = -0.026 ± 0.005, (13) 

which will be used as an independent measurement of A a in the analysis of R^i results in the 
next Sect. 



11 Comparison with previous measurements 

A compilation of measurements of Rm is shown in Fig. 8 (see also Table 6) together with the 
HERA-B results for the three targets. Previous measurements of R^i for atomic mass numbers 
A > 2 were fitted by the function 



Riib' • A , (14) 



with two parameters, where R\m is the value of the fit function when formally extrapolated to 
A = 1 and A a is the power. The E866 value (13) was used as an additional measurement of the 
power in the fit. The fitted parameters are: 

R w = 0.0184 ± 0.0004, A a = -0.030 ± 0.004. (15) 
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Measurements performed with hydrogen and deuterium targets do not follow this simple power 
law dependence, as seen in Fig. 8, where the NA51 measurements [45] for pp and pd interac- 
tions are also displayed. The HERA-B results are in good agreement with the fit of previous 
measurements. 
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Figure 8. A compilation of measurements of R^i as a function of the atomic mass number of the target, A. 
All previous measurements (Table 6) for A > 2 (i.e. excluding NA51) are fitted (dashed line) by the function 
(14) with the E866 value (13) used as an additional measurement of the power in the fit. The HERA-B results 
combined for e and fi are also displayed. 



Table 6. Ratio of ip' and J/tp cross sections in the dilepton channel measured in previous experiments of charmo- 
nium production by protons with momentum, p, interacting with targets with the indicated atomic mass number, 
A. The center-of-mass energy of the interacting nucleons, y/s, is also given. The results of NA50 are listed for 
the "high intensity" samples since they are better suited to study the R^i ratio [33]. The measurements were 
performed in the dimuon channel, with the exception of the ISR experiment (last row) which detected e + e~ pairs. 



pN 


A 


p(GeV/c) 


VS(GeV) 


B'a'/Ba 


Experiment 


PP 


1 


450 


29.1 


0.0157 ±0.0004 ±0.0002 


NA51 [45] 


pd 


2 


450 


29.1 


0.0167 ±0.0004 ±0.00025 


NA51 [45] 


pBe 


9 


450 


29.1 


0.0173 ±0.0004 ±0.0002 


NA50 [33] 


pAl 


27 


450 


29.1 


0.0173 ± 0.0003 ± 0.0004 


NA50 [33] 


pCu 


64 


450 


29.1 


0.0164 ± 0.0002 ± 0.0002 


NA50 [33] 


pAg 


108 


450 


29.1 


0.0157 ±0.0002 ± 0.0002 


NA50 [33] 


pW 


184 


450 


29.1 


0.0153 ±0.0003 ± 0.0002 


NA50 [33] 


P W 


184 


200 


19.4 


0.0180 ±0.0017 


NA38 in [46] 


pU 


238 


200 


19.4 


0.0177 ±0.0022 


NA38 in [46] 


pBe 


9 


400 


27.4 


0.017 ±0.005 


E288 [47] 


pC 


12 


225 


20.6 


0.007 ± 0.004 


E331 [48] 


pC 


12 


225 


20.6 


0.016 ±0.009 


E444 [49] 


pLi 


7 


300 


23.8 


0.0188 ± 0.0026 ± 0.0005 


E705 [32] 


pSi 


28 


800 


38.8 


0.0165 ±0.0020 


E771 [31] 


pAu 


197 


800 


38.8 


0.018 ±0.001 ±0.002 


E789 [34] 


PP 


1 




62.4 


0.019 ±0.007 


ISR [6] 
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The HERA-B measurements of Rj,i for three targets are well fitted by the function (14) with 
the value A a from (15) used as an additional measurement in the fit, resulting in 

R w = 0.0180 ± 0.0006 ± 0.0005, 

A a = -0.029 ± 0.004. ^ ' 

The latter R%-^i value is in agreement with (15). A joint fit of previous (A > 2) and HERA-B 
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Figure 9. A compilation of measurements of R^> in the dilepton decay mode as a function of the center-of-mass 
energy. The measured ratios are rescaled by ^l 029 to compensate the nuclear suppression for targets with differing 
atomic mass number A (see Eq. (16)). The HERA-B results combined for e and fi are shown separately for C, Ti 
and W. The previous measurements from Table 6 (with atomic mass number, A > 2) are presented. The CEM 
expectation [35] is displayed by the dashed line. The NRQCD calculation [50] is shown by the solid line. The 
uncertainty of the calculation is shown by dotted lines. 



measurements, using the value A a from (13) as an additional measurement in the fit, results in 
R w = 0.0183 ± 0.0003 and Aa as given in (16). 

A compilation of measurements of Rm as a function of center-of-mass energy is shown in 
Fig. 9. The measurements were performed with various targets and therefore need to be adjusted 
to compensate for nuclear effects. We assume the ^-dependence given by Eq. (14) with Aa from 
(16) and therefore rescale all measurements by A 0,029 neglecting the error in Aa. The rescaled 
R^i measurements are consistent with a flat energy dependence, in agreement with the CEM [35]. 
NRQCD calculations for R^i [50] show a slow increase with center-off-mass energy. Due to the 
uncertainties of the calculation such a variation of R^> is not excluded, but is disfavored. 



12 Conclusions 

We have performed a study of tp' decays into and e + e~ using a sample of dilepton triggered 

data recorded during the HERA-B 2002-2003 running period. The sample was divided roughly 
equally between the dimuon and dielectron decay modes. The analysis was based on the selection 
of dilepton events with relatively low background contamination and fitting of the dilepton 
invariant mass spectra in the area around the J ftp and ip' signals. The fitted ratio of the numbers 
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of events in the ip' and J/ip peaks was corrected for the J/ip(ip') efficiency, e(e'), to obtain the 
ratio of the ip' and J/ip production cross sections in the dilepton channel, Rm . 

For the analysis, we selected the kinematic domain (2), in which the ip' signal was clearly 
visible. The efficiency ratio, e/e', was mainly determined by geometric factors and, therefore, 
was stable during the running time and reliable for MC calculations. The HERA-B results are 
the first measurement of R^i in the negative xp range. The measured R^i shows little or no ij? 
dependence and is consistent with the previous measurements for positive xf and with CEM and 
NRQCD calculations. The results suggest a slow increase of R^i with increasing px, although a 
flat pt dependence cannot be excluded. 

From the analysis of the angular dependence of R^,', the difference of the polarization pa- 
rameters A for the ip' and the J/ip has been derived: 

AX = A' - A = 0.23 ± 0.17. 

The ratios of ip' and J/ip production cross sections in the dilepton channel combined for the \x 
and e modes and measured for carbon, titanium and tungsten targets are: 

B'a' I Ba (C) = 0.0163 ± 0.0006 ± 0.0005, 

B'a 1 1 Ba (Ti) = 0.0199 ± 0.0024 ± 0.0006, (17) 
B'a' I Bo (W) = 0.0162 ± 0.0010 ± 0.0005, 

respectively. The quoted errors are statistical and systematic. 

Assuming a nuclear dependence in the form of a power law, R-y^i ■ A Aa , and using as a 
constraint the value Act = —0.030 ± 0.004 obtained from previous measurements, a fit to the 
above values yielded: 

R Vljj , = 0.0180 ± 0.0006 ± 0.0005. 

With a joint fit to HERA-B results and previous measurements for different targets (A > 2), 
assuming a power law dependence on the target atomic mass number and using the value Aa 
extracted from E866 results as an additional measurement in the fit, we obtained the following 
values for the fitted parameter: 

R w = 0.0183 ± 0.0003, Aa = -0.029 ± 0.004. 

The HERA-B results (17) indicate no significant energy dependence of R^ with respect to 
previous measurements at lower energies. 

Averaging over all three targets yields the ratio: 

Rp(ji)/Rp{e) = 1.00 ± 0.08 ± 0.04 

This result confirms e — [i universality in dilepton decays of the ip' with an accuracy which is 
better than can be achieved using only current PDG branching ratios for the ip 1 and J /ip dilepton 
decays. The result can be used as an additional constraint to update the PDG branching ratios 
for ip' dilepton decays. 
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